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Abstract 

Matrix model correlators show universal behaviour at short distances. We provide 
a derivation for these universal correlators by inserting probe branes in the underlying 
effective geometry. We generalize these results to study correlators of branes and their 
universal behaviour in the Calabi-Yau crystals, where we find a role for a generalized 
brane insertion. 
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1. Introduction 

Topological strings on Calabi-Yau manifolds provide diverse connections between 
string theory, geometry and random matrix models. In particular it has been shown 
that B-model topological strings on certain non-compact Calabi-Yau backgrounds can be 
described in terms of matrix models [|l|, 0, ||, where the geometry emerges from the 
planar limit of the random matrix integral. Here the matrix model appears as the open 
strings living on compact branes. 

The geometry of these B-models can be effectively analyzed by making use of non- 
compact B-branes [|J. Since these branes have infinite world- volume they are infinitely 
heavy and should be considered as non-dynamical external probes. More specifically, we 
are considering non-compact Calabi-Yau manifolds that are given as the hypersurface 

zw-H(y,x) = 0, (1.1) 

where z,w,y,x G C. The non-compact branes are then parametrized by a fixed point in 
the (x,y) plane that lies on the curve 

C: H(x,y) = 0, 

and they extend in the coordinate z or w. 

The B-model describes the complex structure deformations of the complex curve C 
or, equivalently, of the "Hamiltonian" H(x,y). The variations of the complex structure at 
infinity can be introduced by a chiral boson, <p(x). In a local coordinate patch this chiral 
boson is defined by 

y(x) =d(j)(x), 

and it describes the variation of the curve through its parametrization y(x). 

The target space theory of the B-model on this geometry is the Kodaira-Spencer 
theory [|5j where d(p is the dimensional reduction of the KS field A. The B-branes can be 
thought as defects for the KS field <p(x). The operator ijj(x) that creates or annihilates a 
brane turns out to be a free fermion field and is related to the chiral boson via the familiar 
bosonization formula 

Similarly the field ip*(x) = e~ l< ^' 9a creates/annihilates an anti D-brane. 



In terms of the matrix model, the chiral boson is the collective field of the eigenvalues 
of the matrix. The fermions are basically free fermions, but transform between the different 
patches as wavef unctions with generalized Fourier transformation. We will make use of 
this fact by choosing good coordinates, and later Fourier transforming back to the original 
coordinates. 

We will use the fermionic formulation to study certain correlators in the matrix models. 
In random matrix models, it has been understood for a long time that eigenvalue correlators 
have an interesting behaviour at short scales, called "universality" ||. In particular, the 
joint probablity of n eigenvalues is given by a determinant of a single kernel 

p(Ai, . . . A n ) = det K(Xi, Xj). 

nXn 

In the limit N large, while keeping the rescaled distances N(Xi — Xj) fixed, the kernel K 
takes the form 

where p(X) is the density of eigenvalues at the mean A(Aj + Xj)/2. The formula is called 
universal because it does not depend explicitly on the form of the matrix model potential, 
thus it has the same form for any potential. It has only a functional dependence on the 
potential through the relatively uninteresting scaling factor of the mean density. 

There are many ways to derive this kernel in random matrix models, the usual one 
relies on introducing orthogonal polynomials. With the use of brane insertions, we provide 
a simple derivation by writing the kernel in terms of the correlator of free fermions as 

i((j>(x 1 )-cf>(x2))/gs 

K( Xl ,x 2 ) = (tP(x 1 )r(x 2 )) , (1.3) 

X\ — X2 

where the second equality is by bosonization. Here <p(x) is the chiral boson which is 
related to the geometry associated to the matrix model. In this formula, the use of good 
coordinates, which are single valued, is an essential point. That is, it only takes this simple 
form if the coordinate x parametrizes a unique point on the curve C. It is only in these 
coordinates the branes, or free fermions, can be inserted at a definite position. 

As we will explain, to find the kernel in the usual double valued coordinates, we 
can use Fourier transformation. In fact, we are just applying standard techniques of the 
semiclassical WKB approximation. The kernel will then be given as a weighted sum of 
contributions, coming from the "images" of the brane in the multivalued coordinates. In 



the multicover coordinates, the position of the brane cannot be fixed unambigously, and in 
fact, the brane insertion is best defined via the Fourier transformation. A version of this 
idea first appeared in using methods of conformal field theory. Here we make it more 
explicit and reformulate it in the recent language of brane insertions. 

The universal correlator is an interesting quantity to consider, since it probes the 
geometry at short distances. We can investigate it with the idea of inserting probe branes 
of the geometry. In fact, brane probes can also be used in the more general context of the 
Calabi-Yau crystal §, @, [0. 

We study certain correlators in the crystal which are the 3 dimensional analogues 
of the eigenvalue correlators of matrix models. We find that such correlators can be 
described by fermion insertions in the geometry. These fermions are "generalized" in the 
sense they depend on two parameters. They are the usual one dimensional chiral fermions, 
but inserted at an arbitrary slice of the crystal, giving an additional parameter. Thus the 
fermions are effectively two dimensional objects probing the 3 dimensional structure of the 
crystal. 



The generalized fermions do not correspond to the Lagrangian brane probes in |TC 
They are different objects probing the interior of the crystal. Finding the precise brane 
interpretation is an interesting open problem. 

It is natural to ask whether the free fermion correlators of the crystal also show a 
similar universal behaviour when scaling parameters appropriately. In fact such correlators 



have been computed in the mathematics literature of random partitions [11]] , and a similar 
structure of universal correlators is found. 

When taking a two dimensional slice of the crystal, in the scaling region we find the 
same sine kernel as in the case of random matrices. It is well-known in the mathematics 
of partitions that one can define a suitable probability measure which is the discretized 
version of the measure for Gaussian matrix models. The three dimensional analogue of the 
sine kernel is a more complicated object, written in terms of an incomplete beta-function. 
We conjecture this beta-function is the universal scaling limit of a certain 2-matrix model 



with a unitary measure, which appears in the context of Chern-Simons theory ||12||. and 



the topological vertex [13 
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Fig. 1: The graph of the function y(x) is given by the orbit C determined by 
H(x,y) = 0. In general it has locally several branches: to a single value of x 
correspond the geometric points P\ , Pi , 



2. Branes and the WKB approximation 

To understand the relation between eigenvalues in matrix models and brane insertions 
in topological strings, let us first review some familiar facts about semi-classical quantiza- 
tion. We consider (x,y) as coordinates on a two-dimensional phase space. The curve C 
given by H(x,y) = can be viewed as the orbit that corresponds to the classical ground 
state of the system described by the Hamiltonian H. To such a one-dimensional curve one 
can associate the (global) action variable \x given by the contour integral 

V- = tt f V dx - 
^ Jc 

The Bohr-Sommerfeld rule quantizes this as \i = (N + \)g s , where we identify here and 
subsequently h with the string coupling constant g s . It is useful to think of the Hamiltonian 
H(x, y), and therefore also the curve C, as part of a family of systems parametrized by the 
modulus fx. 

In general, if we want to use x asa coordinate and solve for y as a function of x on 
the curve H(x,y) = 0, there can be several branches yi(x). This is because the line of 
constant x will intersect the curve C in different points Pi, see fig. 1. 

All these points at these different sheets will contribute to the semi-classical wave- 
function ijj(x) associated to the orbit C. The phase of this wavefunction is given by the 
local action 



4>i(x) 



yi(x')dx' 



.I'd 



with xq some arbitrary reference point. The full WKB expression for the wave function 
reads with some overall normalization c 

.(piix) 



^)= c -e(^|) -p 



l- 



9s 



Here all the intersection points Pj contribute. For example, for the harmonic oscillator 
with 

H(x, y) = x 2 + y 2 — Ajj, 

we have two such points. Indeed the wavefunction is well-known to be given in the WKB 
approximation by (for the "allowed zone") 



1p(x) ~ {&JJL — X 2 )' 



4 cos I / a/ Aji — x' 2 dx' — — I 



In the context of matrix models for the topological B-model, the geometric points 
Pi on the curve C correspond to non-compact D-branes on the Calabi-Yau space. The 
variable x appears naturally in the matrix model as the eigenvalue, but will not be a good 
coordinate, since the spectral curve C will always be a multiple cover over the rr-plane 
(a double cover for a single matrix integral). The natural fermion operators i^(x), that 
appear in the matrix model and correspond to creating a single eigenvalue, are therefore 
not represented by a single geometric brane. Instead they are a superposition of branes 
inserted at the inverse images Pi of x, just as we have in the formula for the WKB 
wavefunction. 

In particular for the local geometry with two sheets, where we forget about all the 
^-dependence and therefore can write the local geometry as 

2 2 

y =p , 

we have branes inserted at the points P\ and P2 given by y\^ = ±P- The corresponding 
action functions are <f>i(x) = px and $2 (x) = —px. The wavefunction is naturally given as 
the sum 

If ^ -1/2 \ ■ fP X , n \, . ( PX 7T 

ipix) ~ p ' expz h— +expz 

V \9s 4/ V 9 s 4 

That is, we have in terms of brane insertions the natural identification 



ip(x) = e^i/K-Pi) + e 4 ^(P2)- 
5 
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Fig. 2: Close to a turning point, given by y 2 = x, the good coordinate is y not x. 
Similarly the anti-brane is inserted by 

(Note that this is not the complex conjugated expression.) 

With this translation we can easily compute multi-fermion correlation function. Here 
we should remind us that we will have only contraction of operators that create or anni- 
hilate branes on the same sheet. For example, we have the two-point function 



This can be evaluated in the limit x — > x' as 

sin (y/Ji(x - x')/g s ) 



mx)r(x')) 



J~ji(x — x') 



which gives the famous sine-kernel. 

This description becomes problematic at caustics or turning points where the curve C 
is perpendicular to the x-direction and two of the branches coalesce as in fig. 2. The right 
description at these points is to use the coordinate y instead. The new wavefunction is 
then determined by a Fourier transform 

ip{y) = f dxe lxy ! 9 °iP{x). 



As we will see later this naturally leads to the Airy kernel. 



3. Universal correlator in random matrix models 

We now turn to a detailed investigation of the universal correlators in matrix models 
using the geometric picture of brane insertions just described. Consider an N x N matrix 
model with partition function 

Z= fd&e-h™™ (3.1) 

and let us define the density of n of its eigenvalues as 

p(X 1 ,...Xn)= i -^^(f[Tr5(X k -^). (3.2) 

k=i 

For the random matrix model this describes the probability of measuring n eigenvalues 
at the same time. More precisely, the interesting part of this correlation function is its 
connected part. In the large iV limit the multicorrelator factorizes into one point functions, 
and therefore only its connected component contains new information. 
Introducing orthogonal polynomials 

W(A) 

d\ e 3s P n (X)P m (\) = S nrn 



the joint eigenvalue distribution can be expressed @ , Jl4| as a determinant a single two- 
point kernel 



N n (N-n)l 

AH nxn 



p n (X x , . . . An) = l -, h det K(Xi, A,). (3.3) 



In terms of the wavefunctions 



W(A) 

*n(A.) = P n (A,)e- — 



this kernel can be written as 

K( K A,) = 1 N f ».(A,)».(A,) = E^lg^zlMrl^^MM , (3.4) 

iv '-^ Ai — Aj 

Here we dropped the overall normalization factors, and at the second equality used the 
Darboux-Christoffel formula for the orthogonal polynomials. This formula can also be 
interpreted as a Slater determinant for second quantized fermions 

lf(A < ,A i ) = <* t (A < )*(A i )). 

7 



The eigenvalues of matrix model are usually located along a cut on the real axis. 

When scaling close to a fixed point in the middle of the cut, this kernel satisfies universal 

properties 

sin [NirjK - AXA)] 

K{X - ^ NrriXi-Xj) (3 ' 5) 

in the region where (A; — Aj) ~ 0(1/N). Here A = (Aj + Xj)/2, and p(X) is the single 
eigenvalue density. 

When scaling close to the endpoint of the cut, the kernel has a different expression in 
terms of Airy functions 

Kn ,, Ai(A t )Ai-(A J )-Ai(A J )Ai-(A,) 

K{\i,\j)~ , {S.b) 

Aj — Aj 

where the X's are now suitably scaled variables close to the endpoint. 

Since the form of the kernel does not depend on the specific form of the potential, we 
can restrict attention to the case of Gaussian potential V(x) = x 2 /2. For the Gaussian 
potential, the sine and the Airy kernel can be derived by using the direct asymptotic 
expansion of the wavefunctions, which in this case are Hermite polynomials. We include 
this standard computation in the Appendix. 



4. Universal correlator from geometry 

In the recent context of topological string theory, matrix models are directly related 
to the geometry of the CY. In the limit N — »■ oo, g s — > 0, \x = g s N fixed, the CY geometry 
associated to the matrix model is 

zw + y(x) 2 + {W'{x)f - f(x) = 0. 

Here the quantum deformation f(x) is a polynomial computed from the potential as 

4/j ^ W'{x) - W'(Xi) 



/w = ^E 



N ^ x-X % 

i 

and 

y(x) = d(f>(x) 

is the chiral boson from the collective field 

4>{x) = 2g s Ttlog(x - $) - W'(x) 



The Riemann surface associated with the Gaussian matrix model in this way is given 

by 

H (y, x) = y 2 + x 2 - A41 = 0. (4.1) 

As discussed in |||, the geometry can be probed by inserting non-compact branes. Tak- 
ing an asymptotic patch in the geometry, and considering the worldvolume action of the 
noncompact brane one finds x and y are a symplectic pair of variables with canonical 
commutation relations [x,y] = ig s . In terms of operators, 

d 

ox 

Inserting a non-compact brane in this geometry corresponds to adding a fermion ip(x) = 
e i<j>{x)/9s w ith wavefunction 

*(x) = (iJj(x)) = e ^ cl{x)/9s . (4.2) 

This wavefunction satisfies the Schrodinger equation 

JjU^xJ^^O. (4.3) 

For the Gaussian matrix model these wavefunctions are Hermite polynomials. In fact in the 
Gaussian case these are the same wavefunctions as the orthogonal polynomial wavefunc- 
tions of the previous section. (In the general case the orthogonal polynomial wavefunctions 
do not satisfy the Schrodinger equation, due to more complicated normal ordering issues.) 
In this case, we can simply derive the kernel by finding the wavefunctions from the 
geometry, and inserting in the Darboux-Christoffel formula. Let us first show this for the 
Airy kernel. When close to the branch points, the curve is simply 

y 2 + x = 0, (4.4) 

where we rescaled x for convenience, and put the branch point at x = 0. We can insert a 
brane in a definite position in the y-coordinate which is single valued. In the y coordinate 
the wavefunction solves 

y 2 + tgs-^)V(y) = o, 



so it is given by 

\&(y) = e^ 



l d 



(where we dropped the overall constant.) Transforming back to the x coordinates we get 



*(x) = / dy e l 9s-${y) = g s 1 / 3 Ai 



9s - 2/3 x 



In the large iV limit the Darboux-Christoffel formula becomes 

^(x 1 )^'(x 2 ) - ^(x 2 )^'(x 1 ) 



K(x l7 x 2 ) 



X\ — x 2 



where ^(x) are now the wavefunctions defined from the planar geometry, keeping N large 
and \i = g s N fixed. So from the Darboux-Christoffel formula we immediately have 

Ai(:n)Ai'(:r 2 ) - Ai(a;2)Ai / (a;i) 

K(X 1 ,X2) ~ 



X\ — x 2 



where x\$ are rescaled coordinates measuring distances from the branch points. 

The sine kernel can be derived in a similar way by scaling to a fixed point in the cut 
away from the branch points. In this case the curve reduces to 

y 2 + x 2 — 4/i + 2xx = 0, 

where x is the point we scale to, and x is the distance to that point. Let us work again in 
the y coordinates. If dropping the small x-term, the the wavefunction solves 

'^2 



(y 2 -p(xy)V(y) = 0, 



where p(x) = \fAjj, — x 2 is the classical momentum at x. The wavefunction is a sum of 
5- functions 

^(y) = e~5(y — p) + e~~5(y + p). 

The phases are fixed using the WKB approximation discussed in section 2. Fourier trans- 
forming back 

ty{x) = e 4 e 3s -\- e 4 e s s . 

Choosing the symmetric combination and substituting back in the Slater determinant 
formula we get the sine kernel. 

A more precise way is to keep the small x-term, 

d \ 
y 1 - p(xf + 2ig s x— J *(y) = 0. 

10 



The corrected solution is 



*(s0 =exp 



If 



-p(x) 2 y 



2g s x V 3 
Transforming back to x the solution is again Airy function 



*(x) = (2^ s x) 1/3 Ai 



2^ 

9 2 s 



x 



p(x 



^2 



2x 



(4.5) 



Now we still have to expand in x, which is the small distance to the fixed point x. Doing 
a saddle point expansion of the Airy integral we pick up the contribution from the two 
saddle points 

2/1,2 = ±Vp 2 - 2xx. 

Substituting in the Darboux-Christoffel formula 

sinp(xi — X2) 



K(xi,x 2 ) ~ 



Xi - x 2 



we again get the sine kernel. 



5. Correlator from two brane insertions 

An even shorter way to derive the universal correlators is to insert two branes in the 
geometry. In the single valued coordinates, where the branes can be inserted in definite 
position, the kernel can be expressed as 



K(y x ,y 2 ) = (ip(yi)ip*(y 2 ))- 



(5.1) 



An intuitive derivation of this formula would express the joint correlators in the collective 
field, and then rewrite it with bosonization. However, the x-coordinates here are really 
inconvenient, because they are double valued, and an inconvenient averaging is involved. 
This averaging was first described in 0, and can be understood as summing over the 
images of the brane in the multicover coordinate, as discussed in detail in section 2. 



In the single- valued y-coordinate, we start with the formula ( |5.1| ), which can now 
be interpreted as an insertion of a brane and an antibrane in the geometry. The branes 
inserted are free fermions, consisting of a classical wavefunction, and a quantum part 



11 



So the short-distance correlator in the ^-coordinates is 



e £(.<Kvi)-<Kv*)) 



2/1 -2/2 
Consider now again the geometry close to the branch point, 

y 2 + X = 0. 

The classical wavefunction from the geometry is 

^ {y)=e f I Hy) =e ±y 3 /\ 

Transforming back to the x-coordinates 

J 2/1-2/2 

The Fourier transform is easily evaluated by noting that 

e ix 1 y 1 -ix 2 y2 _ ~ l ( Q _i_ Q \ e ixiyi-ix 2 y2 

Substituting and after a partial integration in rescaled coordinates we arrive at 

WMrM) - Ai ^) Ai/ (^)-Ai(x 2 )Ai'(x^ 

X\ — X2 

Let us expand this formula in the limit where, while close to the endpoint, the distance 
between the two points is very small, so ii = x + x, xi = x — x, where x is the small 
distance. Doing a saddle point expansion of the Airy integral we pick up the contribution 
of the two saddle points y = ±y^. Substituting we arrive at the sine kernel formula 

<V0ni)V*(* 2 )> - ^^- 
x 

This is consistent with the previous expression for the sine kernel, since close to the end- 
point y 2 + x = the density scales as p(x) ~ \/x~. Doing the saddle point expansion of 
the Fourier transform, and picking up the contributions from the two saddle points, we 
automatically introduced an averaging in the double valued ^-coordinates. 
We could also consider the more general geometry 

y m + X = 0. 

In this case the wavefunctions in the y-patch are 

V(y) = exp ' 



g 3 (m + l) 
Fourier transforming in the x-patch this gives 



f i ( y m+1 \ 

y(x) = / dy exp — [ixy + i- — - 

J 9s \ (m + l)J 



The two point correlator and its scaling limits can be similarly derived. 

12 



6. Correlators in crystal melting 

Recently a connection has been discovered between the statistical mechanics picture 
of crystal melting and A-model topological strings on Calabi-Yau manifolds. The crystal 
decribes the toric base of the Calabi-Yau, with lattice spacing of order g s . The temperature 
of the melting crystal is l/g s Q, |J. At very large scales the Calabi-Yau is described in 
terms of classical geometry. Decreasing the distance scale to the string scales, the geometry 
takes the form of a smooth limit shape. At even smaller distances a "gravitational foam" 
picture emerges. 

The structure of the crystal can be probed by putting probe branes in the geometry. 
The effect of adding a single brane probe was first studied in || , where the correction to 
the partition function was discussed. 

Addition of multiple branes and antibranes was considered recently in [TO]. In this 



paper, a statistical interpretation of non-compact probe branes as defects in the Calabi- 
Yau crystal is given, and correlators for inserting branes are derived. These results are in 
accordance to that non-compact brane probes can be thought of as adding free fermions 
with a certain geometric transformation law H, at the B-model side. 

Similarly as for the random matrices, we will investigate certain correlators by in- 
serting branes in the geometry. The correlators we are interested in are the 3D crystal 
analogues of the eigenvalue correlators of random matrices. 

6.1. Partitions on the plane 

Let us first consider partitions on the plane. If we draw the diagram of a partition of 
N, and rescale it with yN, a smooth limit curve of the diagram emerges |f[T |, [15] . This 



is the two-dimensional analogue of the limit shape of the melting crystal. In general, the 
boundary is a graph of a piecewise linear function, consisting from jumps ups and downs. 
Correlators in this context measure the probability of a given pattern of up and down 
jumps. In fact it is enough to consider the probablities to measure the set of down jumps, 
since this already gives the complete information about the shape of the partition. When 
scaling iV large, so that the elements in the set of downs X{ 

■'■'' e [-2, 2] 



such probabilities are given by the determinant of a single discrete sine kernel [|ITJ , [[15 



sm^Xi-Xj) _! / Xi 

K sill (xi,Xj) = f- (pi = cos — = (6.1) 
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Thus the same sine kernel as for random matrices emerges. The important difference 
however is that in this case Xi are discrete, while in the random matrix case we had a 
continuous distribution. Also the scaling factor (j> is different from the random matrix case 
where we had a scaling with the density, but such scalings can always be absorbed in the 
redefinition of the variables. 

The endpoint formula for the random matrix correlator in terms of Airy function also 
has a counterpart in random partitions. 

6.2. 3D Partitions 

We will now study similar correlators in 3D partitions. In the 3D generalization we can 
build a partition from a sequence of diagonal slices, which satisfy the interlacing condition, 

i, !, 0- 



Another way of looking at the partition is viewing it from the top, and projecting 

it to a two-dimensional tiling pattern a. If the partition is located in the positive corner 

of (x,y, z) plane, the position of plaquettes in the tiling pattern is given in the discrete 

coordinates 

t = y-x, 

h = z- -(y + x). 

One can then ask for the probability of measuring a set of fixed plaquettes, at positions 



{(ti, hi) . . . (t n , h n )}. These correlators are given as a determinant of a 3D kernel [16 



KaD((ti,hi),(tj,hj)) = 



dsdw 



(27ri) 2 J| s | =1±e J| w | =lqFe ' s-w s / H+ M+! w; _/ l ._M + i 

Here the contour integral the top/bottom ± signs are valid when t\ > t 2 and and t\ < t 2 
respectively, $^ = l/$3£>, and 

,. / ,>. iim>max(0,-t) i 1 ~ OT I S ) 1 

$3D(«,t) = - rf , 771GZ + -. (6.3) 

This kernel measures the probability of finding two fixed plaquettes in the random tiling 
pattern. In fact, this two-point correlator can be thought of as a non-compact brane and 
anti-brane probe insertion in the geometry. 



The projection is best understood in terms of pictures, see for example [11], [15]. 
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The A-model geometry is related by mirror symmetry to B-model on the mirror 
Calabi-Yau @ 

zw - e~ u - e v + 1 = 



Inserting a single brane, the wavefunction is an eigenfunction of the Hamiltonian 

H(u, v) L(u, q) = 0, 

H(u, v) = q- 1/2 e~ u + e~ 9A - 1, q = e~ 9 ° . 
The one-point function which satisfies this equation is the quantum dilogarithm 



(6.4) 



L(u, g ) = n (l-e"^ n+ ^). (6.5) 

n=0 

Similarly, the one point function for an antibrane is L*(u) = 1/L(u). 
Let us introduce the new variables 



w = e u , 



% h 1 ' ** 
h l = h i -- + —, 



1 \tj 
- + — 

2 2 

In terms of these variables, the kernel can be rewritten as 



hj =h j + ^ + L f 



K ZD {hiX) = 7^-yI [ dudil -^^'^^Kf)^^^ (6-6) 

This is essentially the same form as the two-point free-fermion function of a brane and 
antibrane considered in ||, now Fourier-transformed in the variables (hi,hj). Here <&3D 
and $^" D are the brane and antibrane one-point functions, 

L(u,q) t>Q 

1 i ,1 — 11 _|_ -tn n\ 



1 L(-u 
) L(u- 


+ tg a ,q) 
-tg a ,q) 


1 


-u,q) 



t< 

V ^{~u,q) 

The formula shows that shifting the coordinate £ shifts us to another slice of the crystal, 
parametrized by t. In particular the wavefunction at the slice t = is 

<?> ( r\\ L (u,q) 

$3d(«,0) - 



L(-u,q)' 
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6.3. Brane picture 

It is a natural question what is the precise brane configuration which reproduce the 
3D correlators, and how these branes are inserted in the geometry. To understand this, 
we consider the picture of |l(J where Lagrangian branes were described as defects in the 



crystal. 

Consider now the 3D crystal constructed from a sequence of interlacing 2d partitions, 
ji{t) °. A single 2d partition is a nonincreasing sequence of non-negative integers \x = 
{/xi,//2,- • •}• To build a 3D partitions, such sequences are labeled by an integer t, and 
must satisfy the interlacing condition 

fi(t) < fj,(t + 1) t < 0, 
fjt(t + 1) < fjb(t) t > 0, 

where n(t + 1) > fJ,(t), if 

l*l(t + 1) > fjn(t) > frit + 1) > fMi(t) . . . 

A single two-dimensional slice can be described as a state in the fermionic Fock-space. By 
introducing the set of down jumps of the tableaux and its transpose 

1 

a>i = Hi -i + -, 

t • 1 
bi = Hi -i+ -, 

the fermionic state describing the tableaux consisting of d boxes is given by 

d 



i=l 

Using bosonization, we introduce 

q Lo \ri=q M \v), 
where |/x| denotes the number of boxes, and the creation/annihilation operators 

T ± (z) =e* ±iz) . 



Here we follow the notation of || and |To| . 
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Here <fi±(z) are the positive and negative mode part of the chiral boson <p(z), which is 
related to the complex fermion by bosonization 

ifj(z) =: e*(z) : . 

In this language, introducing a D-brane corresponds to the insertion of a fermion. Dis- 
carding the zero mode part, the D-brane operator can be written as 

y D (z) = TZ 1 (z)T + (z). 

We use here and throughout the paper the standard framing p — 0. 

The partition function of the crystal with the operators can be written as 

oo oo 

Z(q) = (0| J] r +(^) II r-(?- m+ ^) |0>, (6.7) 

n=l m=l 

where q = e~ 9s . With the commutation relation 

r+(z) r_(/) = (i - z/z')- l r- CO r+(*), (6.8) 

it is straightforward to show that the partition function is the McMahon function 

oo 
Z(q) = H(l-q n )- n = M(q). 



71 = 1 



In [[TO]], Lagrangian D-brane insertions were considered. These branes end on the axis, 
and have the geometry 

y = x + a = z + a, a>0, 



where we chose the brane to end at y = a. In |10| it was shown that introducing such a 
Lagrangian brane at a position a = g s (No + 1), corresponds to inserting a fermion D-brane 
operator ^ D (q- ( - N °+^) at the slice t = N + l, 

00 

Z D (q,N o ) = (0\ Hr + (q n -^)x 

n=l 
JVo+1 oo 

H r_(q- m+ ?)y D (q-( No+ ^) l\ r_( ? - m +5) |0). 

m=l m—No+2 
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Commuting through the operators gives the brane one-point function 

oo 

Z D (q,N )=M(q)S{q)]l(l-e-°'( N ° + &q n - h ) 

n=l 

= M(q)S(q)L(g a (N + ±),q), 

where we have expressed the wavefunction in terms of the overall McMahon function M(q) 
and the quantum dilogarithm L(w, q) we introduced before, and u 



oo _. 

^) = n r^- 



q 

These brane insertions depend on the single parameter iVo, expressing the fact that 
the branes correspond to the Lagrangian geometry. Comparing with the brane insertions 
for the 3D correlators in the previous section, we see that our brane insertions are of more 
general kind, since they depend on the two independent parameters t and u. In fact, they 
correspond to the insertion of a D-brane operator &d(z) at an arbitrary slice t. 

Let us insert a fermionic operator ^D(q~^) at the slice t = Nq + 1. Note that N 
and iVo are now independent. This is given by the operator expression 



Z D (q,N,N o ) = (0\ IJr + (<f-2 



2 ) x 

n=l 

N + l 



] [ r_(q- m+ ^ D ( q -( N+ ^) ][ r_( ? " m+ 5) |0). 

m=l m=N +2 



Substituting ^d(z) = T_ 1 (z)T + (z) we obtain 

oo 

Z D {q,N,N ) = (0\ Hr + ( q n -^) x 

n=l 
No+l oo 

II r4 g --+i)r:V (Ar+l) )r + (<r (iV+i) ) J] r_( g — +*) |o). 

m=l m=N +2 

Commuting first the T+(q~( N+ 2') to the right gives a factor 

1 



n a -« 



-N-l\-l 



,_ An+2 L(-^ s (iV+i) + (iVo + l)p s ,g) 



£(g) is a normalization factor with respect to the string answer, so that Z crysta i = £(<?) ^string 

H)|. 
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The leftover TZ (q~^ N+ ^') cancels then the iVth T_ in the sequence. The remainder is 
then the McMahon function, up to an extra factor making up for the missing T_ at the 
place N. This extra factor is 



oo 1 

Y[(l-q n+N )=L(g s (N+-),q). 
Thus 



2 

n=l 



1 ^ W L(-g s (N+±) + (N + l)g s ,q) l J 

Comparing with the wavefunctions appearing in the 3D correlators we obtain 

Z D (q, N, N ) = M(q) ^ 3 d(9s(N + 1), N + 1). (6.10) 

Since we have started with a fermionic insertion at u — g s (N + 1/2) and at the slice 
t = Nq + 1, there is a precise agreement. The branes inserted in the 3D correlator with 
wavefunction $30 (u, t) correspond to the fermionic insertions ^(e 9aU ) at the slice t. Note 
that we have the choice to insert the fermion at a slice t > 0, or at t < 0. Inserting at 
t < precisely reproduces $3.0 (w,£) for t < 0. The branes are generalized branes, in the 
sense that they depend on the two parameters u and t. 

It is interesting to note that the Fourier transform of the wavefunction of the general- 
ized branes $3£>(w, t) is the Wigner function. Let us take t > 0. Using that the generalized 
brane wavefunction is a product of a brane and a shifted antibrane wavefunction, the 
Fourier transform gives 



W(g s t,y) = / duL(u, q)L* (-• u + g s t, q)e 



iuy 



where we recognize the Wigner function. The Wigner function appeared in the recent black 
hole interpretation of the topological string partition function |nj . It would be interesting 
to investigate this connection further. 

6.4- Universality for 3D partitions 

It is interesting to observe that a universal scaling behaviour also appears in the 
case of the correlators of the 3D random partitions. In the same scaling limit as in the 
two-dimensional case, the 3D kernel reduces to the incomplete /3-kernel [JTTf 



K l3 (t,h) = ^- f dz{l 



z y 1 

> z h+t/2+l " 
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where the distances kept fixed are t = ti — tj and h = hi — hj, and r\ is a density scaling 
parameter. This reduces to the plane sine-kernel on two different slices. If we take t = 0, 

or x = y in the original coordinates 

sinfr(/) 
W) = -^-, 

we obtain the sine- kernel in h. Taking h + t/2 + 1 = 0, or x — z — 1, and changing variables 
z->(l-*) 

m ,w 7r(t + 1) 
the kernel reduces to a sine kernel in t' —t+1. Here <p is a redefined scaling variable. This 
means the kernel is likely connected to a two-matrix model. Taking a Gaussian two-matrix 
model with 

"dXdYe-^ x2+l i Y2+v ( x > Y ». 



the universal kernel has the form of K(xi — Xj, yi — Uj) = K(x, y). Reducing to y = or 

x = we have 

sin npx 



K{0,y) 



TVX 

sin npy 



ny 

Thus reducing to the two slices would correspond to restricting to the eigenvalues of only 
X or only Y in the two-matrix model. It is clear that the corresponding two-matrix model 
must have an interaction term, otherwise the full correlator would be just a product of 
the two single-matrix (or single-slice) correlators. This is apparently not the case for the 
,9-kernel. 

The two-matrix model in question is most likely the one related to the Chern-Simons 
partition function [T2\ and the topological vertex |13[ . 



To understand the full kernel, we need the four-point function in this two-matrix 
model, and its scaling limit. Such kernels has been considered in the context of Brownian 
motion |TH] for the usual two-matrix models. Using the context of Brownian motion [ff^ 



could lead to a more precise identification. 

The universal scaling formula (incomplete beta kernel) is a probe of the crystal at 
short distances. Translated to topological strings, it may indicate an interesting short 
distance scaling regime of string amplitudes. 
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Appendix A. Universal kernel from Hermite polynomials 

The universal kernel and its endpoint form can also be computed using orthogonal 
polynomials 

/_V(A) 
dX e 3s P n (A)P m (A) = 5 nm . 

Defining the wavefunctions 

* n (A,) = P n (X l ) e~^t, 

and using the Darboux-Christoffel formula, the kernel (up to a constant) can be expressed 

as 

iV ^—^ Ai — A 2 

For the Gaussian potential we have 

*n(x) = - , H n (-=\e «.. 

(2g s 7r)i2 2y/n\ \^/Zg s J 

where H n are the Hermite polynomials. These are the wavefunctions of the harmonic 
oscillator with h = 2g s . Using H' n (x) = 2nH n _i(x) we can rewrite ( |A.1|) 



as 



"(Ai.AjH 2jy e -. ... ^ — | (A .2) 

where A n is the normalization factor of the Gaussian wavefunctions 

A - 1 



(2g s ir)i2%Vn\ 

We need «(A», Aj) in the limit iV — > oo, gr s — > 0, \x = g s N fixed and large. To find the 
asymptotic expansion of the Hermite polynomials, it is useful to rewrite them in terms of 
parabolic cylinder functions U(a, x) as 

H n (x) = 2 n ' 2 e^ 2 U(-n - \, V2x). 

The asymptotic expansion of the parabolic cylinder function U(a, x) is given in terms of 
the quantity Y = ^/4|o| — x 2 . Since a = —N — 1/2, we expand in the region where a is 

~N 
A* 



large, negative. In this region the argument Xi/ ' yfgl ~ x / — A; ~ viV, so x is moderately 
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large. Using the asymptotic expansion in this regime, and neglecting lower order terms in 
the jj expansion, we find 

U(-N - 1/2, x) = :' , • p= cos - / Y(x')dx'- — 

Substituting in (|A.2j ), and expanding when \ — Xj ~ 0(1/ N) we arrive at 



r^J 



7T(A; - Xj) 

where A = (Aj + Xj)/2, and p(X) is the normalized eigenvalue density 

1 



p(X) = —\4p-X. 
2p% 



The overall normalization factor can be fixed from the definitions, using that /?(A, A) 
p(A). Finally 

w(At, Aj 



sin N7rp(X)(Xi — Aj^ 



iV7r(Ai - A, 



in agreement with (|3.5|) . 



A.l. Endpoint asymptotics 

When scaling close to one of the branch points, A — > ±2y / /I, the kernel has a different 
asymptotic expansion in terms of Airy functions. Let us scale to the positive branch 
point 2^/Jl. At large a the parabolic cylinder functions can be expressed in terms of Airy 
functions as 

j/fexj-a-l-l-rg-i.)^ 

where 

x = 2yf\a\£, t = (4|a|) 2/3 r, 

and 

2 

r=Q(eV? 3 i-cosh- i e) N3 

We need the scaling limit when x — > 2y|a|, that is when £ — > 1. In this limit 

E/"(a,a;) ~ 2-^-5 a r ( - - -o) |a|* Ai(t), 
22 



where t = 2|a|s (£ — 1) mesures the distance from the branch point. Close to the branch 
point then the Hermite polynomials have the asymptotics 

i 
H N (x) ~ 2*e^ \N + ^ 6 T 0^-) Ai(«), 

where 

i 

u = ( A^ + i V (x - V4N + 2) . 
Substituting in (|A.2j ) we find 

vl x Ai(ui)Ai'(u 2 )-Ai(u 2 )Ai / (ui) 

U\ - U 2 

(Here we rescaled u by a \j J~gl as doing the substitution.) 
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